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Abstract 



A new construction is presented for point interactions (PI) and generalised 
point interactions (GPI). The construction is an inverse scattering procedure, us- 
ing integral transforms suggested by the required scattering theory. The usual 
class of PI in 3 dimensions (i.e. the self adjoint extensions of the Laplacian on the 
domain of smooth functions compactly supported away from the origin) is recon- 
structed. In addition a 1-parameter family of GPI models termed resonance point 
interactions (RPI) is constructed, labelled by M. The case M < coincides with 
a special case of a known GPI model; the case M > appears to be new. In both 
cases, the Hilbert space of states must be extended, for M < 0, a larger Hilbert 
space is required, whilst for M > 0, the Hilbert space is extended to a Pontryagin 
space. In the latter case, the space of physical states is identified as a positive def- 
inite invariant subspace. Complete M0ller wave operators are constructed for the 
models considered, using a two space formalism where necessary, which confirm 
that the PI and RPI models exhibit the required scattering theory. The physical 
interpretation of RPI as models for quantum mechanical systems exhibiting zero 
energy resonances is described. 



*E-mail address: C.J.Fewster@amtp.cam.ac.uk 



1 Introduction and Main Ideas 



Point interactions (PI) have long been of interest as solvable models in quantum mechan- 
ics (see [|T3J for an extensive bibliography). Heuristically, they represent Hamiltonians 
with 5-function potentials U H = —A + XS(x)" , although it is well known that such 
Hamiltonians fail to make rigorous sense in dimensions d > 2. Instead, a point interac- 
tion situated at the origin is rigorously defined as one of the self-adjoint extensions of 
the Laplacian on C^°(lR d \{0}), i.e. smooth functions compactly supported away from 
the origin. One of the most useful features of PI is that they represent a leading order 
approximation to the scattering theory of non-point interactions. For example, in 3 
dimensions, the low energy expansion of the S'-wave partial wave shift 5o(k) |I4J] is 



cot S (k) = -— + r k + O(k 3 ), (1.1) 

where L and r are the S'-wave scattering length and effective range respectively, and the 
expansion is valid for spherically symmetric potentials decaying at least exponentially at 
infinity, and for L ^ 0, oo. In d = 3, there is a 1-parameter family of PI, which may be 
labelled {H L | L G HU {oo}}, whose scattering theory is non-trivial only in the .S'-wave, 
where it is given by 

cot<J (A;) = --^, (1.2) 



thus approximating (|1 . 1| ) at leading order. 

The PI described so far suffer from various limitations. Firstly, as noted for example 
by Grossmann and Wu ||, they are restricted only to dimensions d < 3, because —A 
is essentially self-adjoint on C^°(IR d \{0}) in dimensions 4 and higher. Secondly, for 
d = 2,3, they yield non-trivial scattering only in the sector of zero angular momentum. 
Thirdly, they give at best only the leading order approximation to the scattering theory, 
and in the special cases L = 0, oo fail to do even that, for the phase shift for a potential 
with scattering length L = 0, oo does not obey cot5o(k) = 0{k~ l ) at leading order, but 
rather 

cot 5 (k) = pk' 3 + Oik' 1 ) (1.3) 

for L = 0, and 

cot 5 (k) = qk + 0(k 3 ) (1.4) 

for L = oo. To remove these limitations on PI, various generalised point interac- 
tions (GPI) have been proposed, which correspond heuristically to Hamiltonians with 
5-derivative and more general distributional potentials. Models of this type were first 
discussed by Shirokov and were given a mathematical foundation by Pavlov @, |3|, [4j] 
and Shondin || [| (see also @). GPI are not defined on the usual Hilbert space L 2 (1R 3 ), 
but on a larger space: either an extended Hilbert space L 2 (1R 3 ) ©G n Jl|, |5|, Q, |3], or a 
Pontryagin spaceQ n = L 2 (IR 3 ) ©G m 0G m || §, [/J. Due to the presence of negative and 

: A Pontryagin space II is an inner product space which admits a direct orthogonal decomposition 
II = 7i + Q 7i_ into a Hilbert space 7i+ with positive definite inner product, and a finite dimensional 
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zero normed states, Pontryagin spaces do not immediately admit the usual probability 
interpretation of quantum mechanics. However, for GPI models it is often the case that 
there is a positive definite subspace 7i + of II which is invariant under the unitary evo- 
lution generated by the GPI Hamiltonian and therefore has a natural interpretation as 
the space of physical states. 

In conventional treatments of GPI, one decides at the outset how the Hilbert space is 
to be extended, and then constructs the class of GPI which 'live' on this space, either by 
a generalisation of the von Neumann theory of deficiency indices [|J , or by constructing 
separable perturbations of the free Laplacian in accordance with a generalised version 
of Krein's formula || 0, [|- In this paper, we present a new construction of point and 
generalised point interactions, which removes the necessity to determine the underlying 
inner product structure in advance. Rather, the appropriate extension to the initial 
space L 2 (1R 3 ) emerges naturally in the course of the construction. 

Our construction depends on the observation that for a point-like 'potential', the 
scattering phase shifts completely determine the generalised eigenfunctions of the con- 
tinuum spectrum. To see this, note that for a potential compactly supported within 
radius a of the origin, the phase shifts determine the eigenfunctions (up to normalisa- 
tion) for r > a, and therefore if a — > 0, the phase shift in a given angular momentum 
sector determines the continuum eigenfunctions in that sector for all values of the radial 
coordinate.^ Thus, for example, in the S'-wave, we have the radial eigenfunctions Wfc(r) 
at wavenumber k 

/2\ 1/2 

M r ) = (-) sin(£x + So(k)). (1.5) 

We introduce the notation 7i r = L 2 ((0, oo), dr) and 7i k = L 2 ((0, oo), dk) for Hilbert 
spaces of square integrable functions on position and momentum spaces respectively. 
Armed with the generalised eigenfunctions one can define an integral transform 
U : H r — > T~tk so that 

f{r) = (ErVXr) = (-J ^ sin(fcr + 6 (k))f{k)dk. (1.6) 

If U were unitary, then one could immediately define a self-adjoint operator H = U*k 2 U 
for which the Uk{r) would constitute a complete set of generalised eigenfunctions: one 
would thereby have constructed the spectral representation of H. In general, the map- 
ping U is not unitary and so one cannot proceed in this way. For relatively simple 
functions 5o(k), it is, however, possible to determine exactly how U fails to be unitary. 
One can then extend either 7i r or Ti^ (or both) in such a way that U may be extended 
to a unitary mapping U. Thus a self-adjoint operator H = U*k 2 U may be constructed, 
whose generalised continuum eigenfunctions (projected onto the original position Hilbert 
space TC r ) coincide with the Mfc(r). 

Hilbert space TL- with negative definite inner product. The dimension of 7i_ is called the rank of 
indefiniteness. See p8| . 

2 In contrast to the usual treatments of GPI's on Pontryagin spaces, we neglect the possibility of 
distributional contributions at the origin. See the discussion in Section 7. 
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In particular, in the cases studied, it is easy to determine how the Hilbert space 
should be extended. There is a certain amount of ambiguity in our procedure (as in 
other constructions of GPI); however, with the additional requirement of locality, it 
is possible to make natural choices for the various free parameters occuring. We also 
emphasise that our aim is not to construct the most general possible class of GPI, but 
rather to construct at least one GPI with the required scattering behaviour. 

Although certain GPI models have previously been constructed by inverse methods 
|| the analysis usually proceeds immediately from the T-matrix to the discussion of a 
candidate resolvent R(z) whose free parameters are then constrained by the requirement 
that R(z) actually be the resolvent of a self-adjoint operator on the pre-determined 
inner product space. In contrast, the present treatment focusses on integral transforms 
suggested by the scattering data. 

After some preliminaries in Section 2, we illustrate our construction in two cases. 
Firstly, in Section 3, we construct the class of point interactions H L with scatter- 
ing theory given by ( |1.2| ) thereby reconstructing the self-adjoint extensions of —A on 
C£°(IR 3 \{0}), as expected. As is well known, for L > 0, the self-adjoint extension H L 
possesses a single bound state: in the context of our construction, this is manifested in 
the need to extend TCk- For L < 0, there is no bound state, and neither Ti r nor Hk need 
be extended. 

Secondly, in Section 4, we study the class of point interactions on IR 3 modelling 
potentials exhibiting zero energy resonances (infinite scattering length). As noted above, 



the leading order scattering behaviour of such systems ( |1.4j) is not well-approximated 
by any of the usual PI. Instead, we construct generalised point interactions H^ s (M e 
EU {oo}) with scattering theory 

cot 5 (k) = kM (1.7) 

in the S"-wave, and trivial scattering for higher angular momenta. We refer to this class 
of GPI as resonance point interactions (RPI). The special cases M = 0, oo are identified 
with the PI H°°, H° respectively. Apart from this, there are two cases: M > and 
M < 0. For M < 0, we find that H r is extended to L 2 ((0, oo), dr) ©C, although the 
momentum Hilbert space is unchanged. This case represents a subset of the models 
of type B 2 discussed by Shondin [[|. For M > 0, both 7i r and Hk are extended to 
Pontryagin spaces of form L 2 (0, oo) 0(D. This model appears to be new: moreover, it 
is not presently clear how the usual GPI constructions could be used to reproduce this 
model. 

In Section 5, we verify that the PI and RPI models exhibit the required scattering 
theory by explicitly contructing the M0ller wave operators (in a two space setting, where 
necessary). In Section 6, we consider the physical interpretation of RPI models. To do 
this, we employ a general methodology for discussing the 'large scale effects of small 
objects' developed by Kay and the author Hl2| . In particular, we develop a fitting 



formula (analogous to those given in [12]) for matching a given potential V(r) with a 



zero energy resonance to the 'best fit' RPI. This leads to a natural interpretation of the 
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'extra dimension' in the extended Hilbert space of the M < RPI models as representing 
a meta-stable state: part of a wavepacket incident on such an RPI disappears from the 
Hilbert space L 2 (IR 3 ) and the 'missing' probability is stored in the extra dimension 
before being slowly released with a 0(t -1 / 2 ) time dependence. 

It might be objected that neither case L = 0, oo is generic, and that, to all intents 
and purposes, the usual PI suffice to describe the leading order behaviour. However, 
the motivation for the present work arose in a consideration of the scattering of charged 
particles off magnetic flux tubes of small radius ||, in which it was found that the 
scattering lengths for spin-| particles generically take the values or oo in certain angular 
momentum sectors. In consequence, the conventional point interactions (in this case the 
self-adjoint extensions of the operator describing the dynamics in the background of an 
infinitesimally thin wire of flux) fail to describe the leading order scattering theory in 
these sectors. The special nature of this system can be attributed to the fact that it is 
an example of supersymmetric quantum mechanics (see, for example PH). Elsewhere 
TTH, we will construct the appropriate class of RPI for this system. 



2 Preliminaries 

We develop various standard properties of the sine and cosine transforms S and C and 
also demonstrate density of certain subspaces which will be needed in the sequel. S and 
C are defined by 

/2 f 00 / 2 r°° 

- drsmkrf(r) (Cf)(k) = \ - drcoskrf(r) (2.1) 
TC JO V 7T JO 

(the integrals are intended in the sense of 'limit in the mean'). Both are unitary maps 
from 7i r to Tik, with inverses 



(S- l f)(r) 



j 9 poo / 9 roo 

\- dksmkrf(k) (C~ 1 f)(r) = \- dk cos krf(k). (2.2) 

V 7T JO V 7T JO 

Define the operators E± : L 2 (0, oo) — > L 2 (IR) by 

:<o. ^ 

E±(p are the even and odd continuations of ip respectively. For clarity, we also define 

the restriction operator J : L 2 (IR) — > L 2 (0, oo) by (J<p){x) = <p(x) for all x > 0. We 
also define Fourier transformation, T by 

(Ff)(k) = / ^e*7(r). (2.4) 



'2tx 

These operators allow us to express S and C in terms of T: 

(Cf)(k) = (JFE + f)(k) (E + Cf)(k) = (FE + f)(k) 

(Sf)(k) = (JFE^f)(k) (E_Sf)(k) = (FE^f)(k). (2.5) 
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These relationships entail that S and C possess many properties inherited from Fourier 
transformation. The following is a simple corollary of the Paley- Wiener theorem. 

Proposition 2.1 Let f G C °°(0, oo) C H r . Then (Cf)(k) ((Sf)(k)) is an even (odd), 
entire analytic function of k, whose restriction to the real axis decreases faster than 
polynomially at infinity. If g G Cq°(Q, oo) C TCk, then (C _1 (7)(r) ( / (iS~ 1 g)(r) y ) is an even 
(odd), entire analytic function of r , whose restriction to the real axis decreases faster 
than polynomially at infinity. 

Next, we define the subspace V L of L 2 (IR, dk) for L G H\{0} by 

V L = {(1 + (kLff^Tf | / G C °°(IR\{0})} (2.6) 
and the normalised vector ipL G TC by 

Mk)=^y\l + (kL)r i/2 . (2.7) 
We then prove the following: 

Lemma 2.2 V L © {\E + ip L \ \ e&} is dense in L 2 (JR, dk) for any L G 1R\{0}. 

Proof: Suppose <p G L 2 (IR) is orthogonal to T>. Then (1 + (kL) 2 ) 1 ^ 2 ^ is the Fourier 
transform of a distribution supported at the origin. Hence 

for some a a , n. Square integrability then forces to be a multiple of E + ipM- • 
We now define the subspaces X>± of 7^ by 

2?J = {(l + (*L) a ) 1 / a C/|/GCS°(0,oo)} 

2>£ = {(l + (*L) 2 ) 1/a 5/|/GCS o (0 l oo)} (2.9) 



Proposition 2.3 X>+ © {A^l | A G (D} and X>£ are each dense in TCk- 

Proof: Note that C °°(IR\{0}) = £+Cg°(0, oo) © £_CS°(0,oo) and hence that £> L = 
E+V\@EJD L . By Lemma [2^, £> L ffi{A£ + ^ M | A G C} is a dense subspace of L 2 (1R, dfc) 
and admits an orthogonal decomposition into a subspace of even functions and a sub- 
space of odd functions. It is trivial to show that the restrictions of these subspaces to 
1R + must therefore be individually dense in TCk, an d the result follows. ■ 

Finally, the following well known identities are valid for all f(r) G C£°(0, oo) 

k(Sf)(k) = (c^f\ (k); k{Cf){k) = - (s±.fj (k) (2.10) 
For f(k) G C °°(0, oo), we have 

(S-'kf^r) = _A( C -V)(r); (C-'kfXr) = Us^f)(r) (2.11) 
ar ar 
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3 Point Interactions 



In this section, we put into practice the construction sketched in Section 1 by recon- 
structing the familiar class of PI in 3 dimensions. Our starting point is the scattering 
theory described by ( |1.2| ) in the S'-wave for some L e 1R U {oo}, and 6t(k) = for all k 
and all £ > 1. We therefore restrict attention to the S'-wave and attempt to make rig- 
orous the heuristic spectral representation given by (|1.6|), extending the various Hilbert 
spaces involved if necessary. 

In terms of S and C, we define the bounded mapping T L : H r — > 7i k by 

T L = (1 + {kL) 2 )- l ' 2 S - kL{l + {kLf)- l l 2 C (3.1) 
in terms of which equation (|1.6| ) may be re-written as 

f{r) = {n~fW). (3.2) 

In the special cases L = 0,oo, Tl reduces to S and C respectively. Restricting to 
L 7^ 0, oo, we now compute T^Tl- We have 



i c+ c-i L 

l + (kL) 2 \ + {kLf 



~ C \ + (kL)* S ~ S ~\ + {kLY C - (3 ' 3) 

Writing C~ x x +(kV)^ ~ < -*~ ] 1 as an integral kernel, we see that 

2 f°° „ cos kr cos kr' — sin kr sin kr' 1 f°° e lfe ( r +r ) _|_ e -ift(HV) 

dk- 



irJo l + (kL) 2 2irJ-oo 1 + (kL) 2 

= -Le^ r+r '^ (3.4) 
L 

and so we find 

c " TVtkLf C - s "iT(kLf S = \ 1 Xl){xl (3 ' 5) 



where xl is given by 



( 2 V /2 

Xl(t)= — exp(-r/|L|). (3.6) 



Similarly, 



„ 1 kL _ i kL „ sgnL , . . , . 



Hence, we deduce 



-77* t J I_ I Xl)(xl I L > , 
" | I L <0. (3i 
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We now consider TlT£ an d find 

t l ti = I- 



(1 + (kL) 2 ) 1 / 2 
kL 



SC 



-1 



kL 



cs- 1 - 



1 + (A;L) 2 )V2 
1 



(i + [kL) 2 y/ 2 (i + {kL) 2 y/ 2 ' 

On the domain X>_ defined in ( |2.9| ), the last term may be re- written: 
kL „_ , 1 L 



(3.9) 



1 + (ArX) 2 ) 1 ^ (1 + (A;L) 2 )V2 



1 



[l + (kL) 2 ) 1 / 2 dr {l + {kL) 2 y/ 2 
1 _-_a fcL 



1 + (kL) 2 ) 1 / 2 



SC 



[1 + (kL) 2 ) 1 / 2 



(3.10) 



using the identities (|2.10|) and fl2.11|) . Thus we have TlTI = I on T>^, which is dense 
by Proposition 2~3\ Hence 

T L T* L = l. (3.11) 
Including the special cases L = 0, oo, we have thus proved 

Proposition 3.1 For L < or L = oo, Tl is a unitary mapping from TC r to TCk- 

In the case L > 0, Tl fails to be unitary, as it has non-trivial kernel: Tl | xl) = 0. 
Note that this deficiency is restricted to a 1-dimensional subspace. We can compensate 
for the 'missing probability' by extending TCk to TLk ©C (with inner product (/ © a \ 
g © (3) = (f | g) ip- + oc/3) and then defining the mapping Tl by 



Tl = [T L , (Xl |] : H r 
f 



Hk ©0 

T L f © (Xl I /) 



(3.12) 



We then have 



Proposition 3.2 For L > 0, Tl is a unitary mapping from H, r to Tik ©0, with inverse 



-i 



T L (f®a)=Tlf + a\ X L). 



(3.13) 



Proof: Note that Tl is a surjection onto Hk as a consequence of fl3.11| ). The result then 
follows immediately from the definition of Tl and (|3.8| ). ■ 

Our extension to the momentum Hilbert space will, of course, carry the interpretation 
of a bound state with normalised eigenfunction xl(t) E H r - Note that the mapping Tl 
is just one of a 1-parameter family of physically equivalent possible mappings T[ = 
Tl, e l5 (XL | corresponding to re-phasing \l- 
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We may now proceed to define the point interaction Hamiltonians. In the case L < 0, 
and also in the special cases L = 0, oo we have seen that the mapping JF L is unitary, and 
so we may immediately define a self adjoint operator by 

h L = T* L k 2 T L (3.14) 

with domain D{h L ) = TIDik 2 ). In the case L > 0, however we must use the unitary 
mapping Tl, defining 

h L = f L ~\k 2 ® E L )f L (3.15) 

with domain D{h L ) = jF~£ l (D{k 2 ) ©(D), where El G IR is arbitrary. We have restricted 
our operator to be diagonal in the momentum representation in order to ensure that 
the continuum eigenfunctions are still given by the Uk(r). Note that our construction 
does not uniquely determine the energy El of the bound state - any El results in an 
unbounded self-adjoint operator whose continuum eigenfunctions have the required form. 
In order to remove this ambiguity, we impose the additional requirement of locality, in 
the form of a requirement that h L should agree with —d 2 /dr 2 on C^°(0, oo). Clearly this 
will automatically restrict us to self-adjoint extensions of —d 2 /dr 2 on this domain. 

Proposition 3.3 El = — \L\~ 2 is the unique value for which h L is local. 



Proof: From Proposition [2J], it follows that JF L C^°(0, oo) C D{k 2 ) and hence C£°(0, oo) C 



D{h L ). For / G C^°(0, oo), we compute 

h L f = nk 2 FLf + E L \ X L)(XL\f) 

= -^Lf" + E L \ X L)(XL\f) 

= -f"+\XL)(EL(XL\f)-(XL\f")) 

= -f" + (E L + \L\- 2 )\xL)(XL\f) (3.16) 

where we have used the fact that k 2 TLf = —Fhf" for / G C^°(0, oo) and also that 
(— d 2 /dr 2 \ C g°(o,oo))*XL = — \L\~ 2 Xl- The result follows immediately. ■ 

We now determine the domain of h L explicitly. This is, of course, well known; we 
discuss it here only to show how it may be derived within the terms of our construction. 



Theorem 3.4 h has domain 

D[h L ) = {<p\<p,<p' E AC loc (0, oo); <p" G L 2 (0, oo); <p(0) + Lcp'(0) = 0} (3.17) 

where the boundary condition y?(0) + L(p'(0) = is to be interpreted as <p(0) = for 
L = 0, and <p'(0) = for L = oo. 

Proof: We give details for the case L < and indicate how the proof is modified for the 
remaining cases. Note that C^°(0, oo) C Tlk is a core for k 2 , and hence V = JF^C^O, oo) 
is a core for h L . Any / G T> may be written 

/ = [S- 1 -C- x kL)^ (3.18) 
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where <p — (1 + (fcL) 2 ) _1 / 2 jF L / e C^°(0, oo). By Proposition ^J] and identities of the 
sine and cosine transforms, we find that f(r) = g(r) — Lg'(r), where g(r) = S _1 ip is 
odd, entire analytic, and decreasing faster than polynomially as r — > oo. In particular, 
/(0) + Lf'(0) = for all / G T>. Moreover, there is at least one fo € T> for which 
/ (0) ^ 0, for otherwise we would have (C -1 Jfe(l + (kL) 2 )- 1 ' 2 ^) = for all r\ in the 
dense set Cq°(0, oo), obtaining a contradiction. Now h L agrees with —d 2 /dr 2 (acting in 
the sense of distributions) on V, and is symmetric on this domain. Thus 

D(h L ) = D{{h L \ v )*) ={<p\(<p\ f") = (<P" I />, for all feV} (3.19) 

and the required result follows easily. The cases L = 0, oo may be treated similarly, 
using l S~ 1 C^°(0, oo) and C _1 C^°(0, oo) as cores. For L > 0, the appropriate core is 
P = JP£C °°(0,oo)©{AxL(r) | AeC}. ■ 

We can assemble the full PI Hamiltonian acting on L 2 (1R 3 ): with respect to the 
decomposition 

oo 

L 2 (IR 3 ) = 0L 2 (IR + ,r 2 dr)®/Q (3.20) 

1=0 

where /Q is the subspace of L 2 (S 2 , dfl) spanned by Yg-i, ■ ■ ■ , H L is defined by 

oo 

H L — U*h L U ® I © U*hiU ® I (3.21) 

£=1 

where U : L 2 (1R + , dr) — > L 2 (IR + , r 2 <ir) is the unitary operator (Uf)(r) = r/(r), and ^ 
(£ > 1) is the unique self-adjoint extension of -d 2 /dr 2 + £(£ + l)/r 2 on C^°(0, oo) C 7^ r . 

To summarise, we have seen how the usual class of point interactions may be con- 
structed from a consideration of the required scattering theory. 



4 Resonance Point Interactions 

We now construct the class of RPI models by a similar method to that used in the 
previous section. In this case, the relevant heuristic spectral representation is again given 
by (|1.6| ), but with 5o(k) now specified by Ql.7|) for some Me RU {oo}. Accordingly, we 
consider the bounded mapping Tm defined by 

T M = (1 + (kM) 2 )~ 1/2 C + kM(l + (kMf)- 1,2 S. (4.1) 

In the special cases M = 0, oo, Tm reduces to C and S respectively. Hence in these 
cases, Tm is unitary from H r to Hk and the operators h^ s = T M k 2 TM are well defined 
self-adjoint operators. Comparing with the previous section, we see that h® es = h°°, and 
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We now consider the remaining cases. Firstly, we compute Tj^Tm- We have 

t^Tm = l + c-\ + l [kM)2 c-s- 1 i ^ kM)2 s 

n i kM ^ ^ i kM /A , 

+c TT(W 5 + 5 TTPF C (4 ' 2) 

and so, by the arguments used in the previous section, 



q-* t-r j 1+ I Xm) (xm I M > , . 



which should be compared with equation ( |3.8|) . 
Computing TmT m , we have 

1 1 /cM 



1 + (A;M) 2 )V2 (i + (kM) 2 ) 1 / 2 

kM , 1 . . 

, 71~77Tovi79 (4-4) 



(i + (kMfyi 2 (i + (kM) 2 y/ 2 

On the domain defined in Q2.9|) , the last term may be re- written: 

7<SC 7^ 77T7T777 = TT^T^T^C — C' 1 



2 



(l + (kM) 2 y/ 2 (l + (kM) 2 y/ 2 (l + (kM) 2 ) 1 ' 2 dr (l + (kM) 2 y/ 

~ M * „ (4-5) 



(i + (kMyy/ 2 (i + (kM) 2 y/ 2 

Proposition pT3| shows that T> + © {A-^Af | A G C} is dense in Hk, so we require the action 
of T M T M on ^m- Explicit computation using the standard results 



yields 



y/2ip M (k) M > 



and 

(T^M)(r) = -^(1 + sgnM) XM (r). (4i 



from which we obtain 



TmT m i/)m = < M < o. ^ 4 - 9 ) 



Assembling our results, we have 



r i+iM&M i m>o 

r ^-| I-I^X^I M<0. (410) 
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We now modify Tm in order to obtain a unitary operator. First note that the M < 
case is analagous to the L > case for JF L , except that here, it is T M T M rather than 
T* l Tl which fails to be the identity. Accordingly, we extend the position Hilbert space 
to TC r ©C (with the obvious inner product) and define the mapping Tm by 

T\i : TL, ©C — > Hi, 

f®a — > T M f + a\^ M ). (4.11) 
T M is unique up to re-phasing of ipM- We then have 

Proposition 4.1 For M < 0, Tm is a unitary mapping from 7Y r ©C to Tik with inverse 

Tuf = [Tmi (V>m I] / = © (V>m I /)• (4.12) 
We may therefore construct a self-adjoint operator h^ s on 7i r ©C by 

C = t M ^ 2 t A/ (4.13) 

with domain D(h^ s ) = T^ I 1 D(k 2 ). We will return to the physical interpretation of the 
'extra dimension' in Section 6. 

We now examine the properties of h^ s . By the same arguments as in the case of L > 
PI, the domain of h^ s includes the (non-dense) subspace V = {(f) © | G C£°(0, oo)}. 
Moreover, one can easily see that 

C(0©O) = T M k 2 TM<f> © (4>m | k 2 TM<p) 
= -T M T M (t>" © (Vm I - W> 

= -</>"© (4.14) 

for G C^°(0, oo) where we have used the fact that T M ipM = for M < 0. Hence, h^ s 
is a self-adjoint extension of the non-densely defined operator 

d 2 

- — ©0 onV C L 2 ((0,oo),dr)©(D (4.15) 

The /i£f s are therefore local, and belong to the class of models considered by Pavlov || 
and also by Shondin In the nomenclature of 0, the /i^£ s for M < form a subset of 
models of 'type B 2 : . 

To describe the domain of hf* s explicitly, we use a similar argument to that employed 
in the PI case, identifying V = 7^ C£°(0, oo) as a core for h^ s . If g(k) G C£°(0, oo), 
then (ip M I = \M\^ 2 (T M g)(0), and so £> may be written 

V = {/ © IMI 1 / 2 /^) | / G T^CHO, oo)}. (4.16) 

Any function / G T M C£°(0, oo) may be written /(r) = g(r) — Mg'(r), where #(r) is 
analytic, even, and decreasing faster than polynomially as r — > 0. As a consequence, 
f(0) = -M/"(0). On £>, has action 

C/©|M| 1/2 /(0) = -f»@-\M\ x l 2 f"{<S) 

= -f"Q)-\M\~^ 2 f\0) (4.17) 
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The domain of h^ s is equal to D((h^ s \x>)*)- We have 

I *&(/e |M|VV(0))) = (-_^|/)+Ro)f(o)-7(o)/(o) 

-^ilMl- 1 / 2 /'^) (4.18) 



for any such that cp" e L 2 (0, oo) in the sense of distributions. Note that —cP/dr 2 is 
no£ symmetric on T^C^°(0, oo). Hence, if y?©$i G D(h^ s ), with ^^©$1 = -<//'©$ 2 , 
then 

7(0)/(0)-R0)/'(0) = $ 1 |M|- 1 / 2 f(0) + $ 2 |M| 1 / 2 /(0). (4.19) 



for all / G 7^C£°(0, 00). Moreover, /(0) and /'(0) are independent on this domain. The 
following is then immediate. 

Theorem 4.2 The domain of h^ s is 

D(h% B ) = {tp®$\<p,<p' e Aa o e(0,oo); </gL 2 (0,oo); $ = IMI 1 / 2 ^)} (4.20) 
u>it/i action 

h&bp © IMf/Vtf))) = V © -\M\-^'(0). (4.21) 



In contrast to the usual PI case, the boundary condition corresponding to ( 4.20| ) is 
energy dependent: the eigenfunction equation at energy k 2 is equivalent to the equation 
— if" = k 2 (p, with boundary condition 

k 2 M<p(0) = <p'(0) (4.22) 

and the generalised eigenfunction at energy k 2 is therefore 

2 . 1/2 



//,, = I -J sin(fcr + S (k)) © |M| 1/2 sm5 (k) (4.23) 

with 60 (k) given by fll.Tp . 

As a result of our construction, the spectral properties of h^ s are easily identified: 
h^ s has purely absolutely continuous spectrum cr(h^ s ) = cr ac (h^ s ) = [0, 00). 

Turning to the case M > 0, we note that our original choice of continuum eigenfunc- 
tions form an 'over complete' set. We can remedy this by extending to a Pontryagin 
space, the heuristic motivation being that this will allow us to subtract off the 'excess 
probability'. It is only possible to construct unitary mappings between Pontryagin spaces 
with the same rank of indefiniteness, so we must extend both position and momentum 
spaces to Pontryagin spaces. As before, we find that the failure of unitarity is located in 
a one dimensional subspace, which suggests that we choose Pontryagin spaces of form 
II = L 2 (0, 00) 0(D, where the indefinite inner product is given by 

(fea\geP)u = (f\9)»-&P. (4-24) 
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We define the mapping T M , unique up to re-phasing of ipM and xm, by 

f M :H r e& — ► H k G& 

fea — > (T M f -a\^ M ))e(-(x M \ f} L *+aV2). (4.25) 

A short computation proves 

Proposition 4.3 For M > 0, T M is a unitary mapping from H r 0© to Hk 0©, with 
inverse 

f M l :H k Q€ — ► H r Q(D 

fea — > (T^f + a\xM))e(^M\f)L* + aV2) (4.26) 

Hence, we can construct the RPI Hamiltonian for M > on L 2 ((0, oo), fir) 0©: 

C = ^m 1 (* 2 e £m) Tm (4.27) 

(We have again chosen the Hamiltonian to be diagonalised by T M .) As in the case of 
PI for L > 0, a single free real parameter Em is introduced by our procedure. Em is 
interpreted as the energy of a negative-normed eigenstate | xm) O v^2- Such a state is 
undesirable in a quantum theory. Accordingly, we decompose the position Pontryagin 
space n r into the orthogonal subspaces TC + and £ (which are independent of the choice 
of Em) given by 

n r = n+ec 
n + = f M \H k ev) 

C = T M 1 (0eC) = {A(|xM)ev / 2)|AGC} (4.28) 

where we use the notation in contrast to the decomposition n r = TC r 0©. 

Due to the diagonal structure of hf^ s in the momentum representation, hf^ s respects 
this decomposition 

C = h^ + QE M . (4.29) 

Moreover, 7i + is an intrinsically complete, positive definite subspace of n r of unit co- 
dimension, orthogonal to the negative-normed eigenstate found above; it may therefore 
naturally be identified as the space of physical states. The operator + and the space 
7i + are independent of the value of Em, so to some extent this parameter has no physical 
meaning. However, if one wishes the operator h^ s to be local on the Pontryagin space 
n r , one can identify the value Em = — |M|~ 2 as the unique value compatible with this 
requirement by similar arguments to those used in the PI case for L > 0. 

We now describe the domain of httL , in more detail. We have D(h^ + ) =f M ' 1 (D(k 2 )G 
0) C H + . The domain V = 7^ 1 (C^°(0, oo) 0) is a core for h^ + , and may be written 

V = {/ \M\^ 2 f(0) | / G T M C™(0, oo)}. (4.30) 
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As before, we note that any function / e TJ^Cq°(0,oo) may be written f(r) = g{r) — 
Mg'(r), where g(r) is analytic, even, and decreasing faster than polynomially as r — > 0, 
and that /'(0) = —Mf"(0) in consequence. The action of h™ s + on V is given by 

<+(/ e \M\ l ' 2 fm = -/" e \M\-^f(o) (4.3i) 

An exact analogue of the argument used in the M < case then yields 

Theorem 4.4 The domain of h™ s + is given by 

D{h^ + ) = {<pe$\(p,<p' e AC loc (0,oo); <p" e L 2 (0,oo); $ = |M|VV(0)} (4.32) 
with action 

O 8 = V 6 IMI-Vy^). (4.33) 

We therefore have the full domain D(0 = D(hf cs + )QC 

DK S ) = {( V 9e$)+a(|xM)ev / 2) |^,^e AC loc (0,oo); 

G L 2 (0, oo); $ = \M\ 1/2 ip(0); a E <D} (4.34) 

and action 

C(v9 e $) + a(\ X m) e V2) = e iMr^V'Co) + «s M (| X m) e v^). (4.35) 

As before, the /i£f s exhibit energy dependent boundary conditions: solving the eigen- 
value equation at energy k 2 yields solutions obeying 

k 2 Mip(0) = <p'(0). (4.36) 

The spectral properties of h^ s are as follows: cr(h^ s ) = cr pp (h^ s ) U <7 ac (/^£,), with 
ffp P (C) = {Em} and a^hfj = [0,oo). 

We can assemble the full RPI models on R 3 as before: defining H^ s by 

oo 

H™ s = U*h™ s U © I © U*h t U © I (4.37) 

e=i 

where, for M < 0, U(f © a) = r/(r) © a, and for M > 0, [7(/ © a) = r/(r) © a. 

5 Scattering Theory 

The original aim of our construction was to produce PI and GPI Hamiltonians with a 
given S'-wave phase shift. It is therefore expedient to check that the models described 
above actually exhibit the required behaviour. The usual method of demonstrating the 
existence and completeness of M0ller wave operators for PI and GPI models is to show 
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that their resolvents are trace class perturbations of that for the free Laplacian, and 
then to apply Kuroda-Birman theory ||16|| . Instead, we will use a method which builds 
on our construction and yields the M0ller operators (and their completeness) directly. 

We work in the S"-wave only, and employ a two space setting: let B be self-adjoint 
on TCi, A be self-adjoint on 7i 2 and J be a bounded operator from Hi to 7i 2 . Then the 
M0ller operators Q^(A, B\ J) are defined by 

^,5: J) = lim e iAt Je~ lBt P ac (B) (5.1) 

t — >=foo 

and are said to be complete if the closure of Ranf2 ± (A, B; J) is equal to RanP^A 
Our results in this section follow from 

Proposition 5.1 U_ t CS~ l U t — > ±zl as t — > =Foo, where U t is multiplication by e _jfc2 * 
on Tik- 

Proof: For any u{k) G Cq°(0, oo), we compute 

\\U- t CS~ l U t u(k) =f iu(k)\\ 2 = ((((T 1 ± iS- l )U t u{k)\\ 2 



2 r°° 



dr 



dke i{±kr - k2t) u(k) 



2 



(5.2) 



7T JO 

and note that the last expression vanishes if t — > =Foo by (non)-stationary phase argu- 
ments (see the Corollary to Theorem XI. 14 in ||16[). ■ 

Using this result, we see that if Q = cos5o(k)S + sm5o(k)C, then 

U-tQS^Ut -> e ±l5 °^ (5.3) 

as t — > =poo. It is then easy to construct the M0ller wave operators for the PI and RPI 
models. For the PI case, we take H\ = 7i 2 and J to be the identity, writing Q ± (A,B) 
for the wave operators. 

Theorem 5.2 Q ± (h L ,h°) exist, are complete, and given by 

fL±(h L , h°) = F* L e ±i5o{k) S (5.4) 

where S (k) is given by ( \T!\ ). 



Proof: For L < 0, L = oo, the existence and form of the M0ller operators is immediate 
from the above, and the definition of h L as JF£/c 2 jF L . Completeness holds because all 
three factors in ( |5.4| ) are unitary, and hence Ranf2 ± (/i i , h°) = 7i r = RanP ac /i L . For 
L > 0, we have 

e ihH e -ihH = jr lU t jr LS -l UtS + e -m~H | ^ ^ | £-1^ (5 5) 

The second term vanishes as \t\ —>■ oo by another non-stationary phase argument, and 
the required result follows because RanjF£ = RanP ac /i i . ■ 
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Theorem 5.3 For M < 0, Q ± (h^ s , h°; J) exist, are complete, and given by 

n ± (h^,h°)=f M 1 e ± ^S (5.6) 
where S (k) is given by and J : 7i r — > 7i r © (D is defined by Jf = f © 0. 



Proof: The existence and form of the operators is due to Proposition [5J] and the following 
observation along with the form of T M and the fact that T M J = T M © 0. Completeness 
holds because RanT^ 1 = Ti r = RanP ac /^f s . ■ 

For the case M > 0, our two spaces are TC r and H+, the physical state space. 

Theorem 5.4 For M > 0, Q ± (h^ s , h°; J) exist, are complete, and given by 

tt+ih™, h°) = t^e^WS (5.7) 

where So{k) is given by flOU , and J : H r —> H r ©C is defined by Jf = P+(f © 0), wit/i 
P + the orthogonal projector onto 71+ . 

Proof: The argument proceeds as before, completeness holding because the wave oper- 
ators are isometries from TC r to 7i + = P ac h^ s . ■ 

We conclude that our construction does indeed yield the required scattering theory, 
and also (because of the way in which the PI and RPI Hamiltonians were defined) that 
complete M0ller operators may easily and explicitly be determined. 



6 Physical Interpretation 



In this section, we show how RPI models may be used to model Schrodinger opera- 
tors H = —A + V, where V is smooth, spherically symmetric, compactly supported 
within radius a of the origin, and exhibits a zero energy resonance. Our methodology is 
analagous to that developed in 
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in which the non-resonant case is discussed. Here, 
we develop a formalism (the fitting formula) for selecting the 'best fit' RPI for such 
operators. The range of energies for which the approximation is valid can be determined 
by a 'believability' analysis analagous to that developed in |12|. We will not do this 
here. 



The S"-wave wavefunction at wavenumber k obeys 




2 + V(r) > u k (r) = k 2 u k (r), 



(6.1) 



where u(r) obeys regular boundary conditions, i.e. u(0) 
within radius a of the origin, the phase shift is given by 



0. Because V is supported 



cot So(k) 



ka sin ka + D(k) cos ka 
ka cos ka — D(k) sin ka ' 



(6.2) 
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where D(k) = au' k {a) /u]~(a). Expanding D(k) in powers of (ka) 2 , D(k) = D +Di(ka) 2 + 
0((A;a) 4 ), and substituting in (|6.2|), we find 

,, (h , Do + (1 + Di-Do)(ka) 2 + O((ka)*) 
COt 6 ° {k) = (l-D )ka + O((kar) ■ (6 - 3) 

This then leads to the expression L = a(Do — 1)/Dq for the scattering length. Clearly, 
a zero energy resonance (infinite scattering length) occurs when D = 0. In this case, 
the leading order behaviour of 6o(k) is given by 

cot 5 (k) = (1 + Dx)ka + 0((ka) 3 ). (6.4) 

Comparing with ( |1.7| ), we find that the leading order approximation to the dynamics is 
given by the RPI H^, with M = a(l + Di). Thus it suffices to compute D\ for the 
potential of interest. To do this, we use the well-known formula (see, e.g. |i~7| , p. 236) 

Dik x ) - D(k 2 ) = -((M 2 - {k 2 af) S>k ^ k ^ T (6.5) 

au kl {a)u k2 {a) 

Using this to expand D(k) about k = 0, we find 

= _a-g«.(r)'Jr 
u {a) 2 

and thus arrive at the fitting formula (cf. (12[ 

V Mo(a) 2 / 

The best fit RPI can therefore be computed in terms of the zero energy solution to fl6.1| ). 
In addition, the labelling parameter obeys the bound 

-oo<M<a. (6.8) 

Moreover, this bound is best possible: for any M in the above range, one can clearly find 
a smooth function Uo(r) satisfying regular boundary conditions at the origin, uq constant 
for r > a and such that ( |6.7| ) holds. Then the potential defined by V(r) = u'l(r)/uo{r) 
has infinite scattering length, and scattering theory approximated at leading order by 

ttM 

res' 

It is interesting to note that the RPI models themselves obey ( |6.7|) in the following 
sense: the zero energy generalised eigenf unction for M < is given by 

m = 1©|M| 1 / 2 . (6.9) 

Interpreting the integral in ( |6.7| ) for elements / © a of H r ©C as J a \ f{r)\ 2 dr + |a| 2 , the 
right hand side of equation (|6.7|) is then equal to M. Similarly, for M > 0, we interpret 
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the integral for / a G 7i r 0(D as J Q a \f(r) \ 2 dr — \a\ 2 and use the generalised zero energy 
eigenf unction 

Mo = ie|M| 1/2 (6.10) 

to yield the value M. 

We briefly consider the interpretation of the extension to the Hilbert space 7i r in 
the case M < 0. From the fitting formula, it is clear that M is negative if and only 
if the mean square value of Uo(r) within r < a exceeds Uq{o) 2 . Hence this case is 
used to model Schrodinger operators whose zero energy generalised eigenfunctions are 
peaked inside the interaction region. This is characteristic of resonant behaviour and 
corresponds to a physical picture of a particle being detained inside the region, before 
being gradually released. If the radius of support is shrunk to a point, the particle must 
be completely removed from the space in order to model this process. It is therefore 
natural that the Hilbert space be extended in this case. 

One may study the time behaviour of the proportion of a wavepacket in the extra 
dimension. Starting with a normalised vector © 1, the evolved packet is given by 
'T-^ 1 e~ %k t ipM- The proportion remaining in the initial state at time t is 

f(t) = (i>M\e- ikH i> M ) 
2 | M | r 



Ti Jo 1 + (kM) 
which may be approximated using the method of stationary phase to give 

<p{t) ~ ^M e --/4 + 0(r i) ( 6 12 ) 



Tit 

as t — > oo. It is interesting that the decay is not exponential. The M > case does 
not admit such a simple interpretation due to the redefinition of the space of physical 
states. 

As a simple example of the use of the fitting formula, we consider the example of a 
square well with a zero energy resonance. Setting V(r) = — {{n+^)n / a) 2 (n = 0, 1, 2, . . .) 
for r < a and V(r) = for r > a, we find Uo(r) = sin(n + ^)nr/a. Equation ( |6.7| ) then 
yields M = a/2 for any value of n. 



7 Conclusion 

We have seen how PI and GPI models may be constructed and studied using integral 
transforms suggested by the scattering data. In addition, we have constructed a new class 
of RPI models (the case M > 0). The RPI models obey energy dependent boundary 
conditions and are defined on a Hilbert or Pontryagin space which extends the usual 
space of states. In the Pontryagin case a positive definite invariant subspace may be 
constructed which may be interpreted as the physical state space. 
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It is instructive to compare our methods with those usual in constructions of GPI 
models on Pontryagin spaces 0. One starts with the desire to include wavefunctions 
with distributional terms in the space of states. Accordingly, choosing a dense subspace 
T> C Tir and an anti-linear functional u (more generally m such functionals) on T> 
{uj 7i r ) one constructs a linear space Ti generated by elements 

i) = i) + \uo (7.1) 

with tj)Q E T> , A G C The next step is to extend the Hilbert space inner product (• | •) 
from P to a sesquilinear form (•, •) on Ti by defining 

{^uj)=uj^ ) tp eV (7.2) 

(recall that uj is anti-linear). This leaves (uj,cu) as a free parameter which is usually 
fixed on the basis of a physically inspired renormalisation of certain divergent integrals. 
The form (•, •) is in general indefinite; however, there is a closely related positive definite 
inner product which is used to form the completion TiofTi. With respect to (•,•), Ti is 
a Pontryagin space Ti r ©C ©C. One can then seek the GPI models which 'live' on Ti. 

The present treatment has various advantages over this procedure. Firstly, one need 
not determine the distribution uj in advance; moreover the form of the Pontryagin space 
is suggested naturally by the over completeness of the generalised eigenf unctions ( |1.5|) . 
Indeed, distributions need never explicitly occur in our construction; although it is clear 
(e.g. from the domain of definition of h^ s ) that the component of the wavefunction in 
the 'extra dimension' carries a distributional interpretation. Finally, it is not clear that 
the usual procedure can encompass the M > RPI models without modification, as the 
natural Pontryagin space is Ti r CD rather than Ti r © C C. It is possible that these 
models are restrictions of GPI models in the larger space; however, it seems more likely 
to us that a better starting point would be to choose a non-dense test space V consisting 
of vectors orthogonal to a 1-dimensional subspace. It is not presently clear exactly how 
this would be implemented, nor what the appropriate distribution uj would be. 

The construction presented here has so far only been employed for interactions whose 
scattering is non-trivial only in the S'-wave. One of the great strengths of the usual 
constructions is that the generalisation to higher angular momenta is relatively straight- 
forward. It would be interesting to extend the present treatment to this case. 

Finally, it is of interest to understand whether RPI models can arise as limits of 
sequences of Schrodinger operators with potentials of compact support decreasing to 
the origin, e.g. in the spirit of [12] in which the usual class of PI is exhibited as strong 



resolvent limits of such sequences (see |L3| for a treatment using sequences of scaled 
potentials in the norm resolvent topology). We hope to address these issues elsewhere. 
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